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Overview 


■  Purpose:  Convey  how  weather  forecasting  models 
incorporate  diverse  observational  data 

■  Governing  Equations 

■  Introduction  to  Meteorologicai  Observations 

■  Introduction  to  Numerical  Weather  Prediction  (NWP) 

■  Methods  of  Data  Assimilation 

■  Optimal  Interpolation  (01) 

■  3-Dimensional  Variational  Assimilation  (3DVAR) 

■  4-Dimensional  Variational  Assimilation  (4DVAR) 

■  Kalman  Filter 
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Atmospheric  Governing  Equations 


Conservation  of  Momentum 
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Atmospheric  Governing  Equations 


■  Conservation  of  Momentum 

—  =  -i-V/7  +  g  +  F^-2ilxU 
dt  p 

Geostrophic  Balance:  Coriolis  acceleration  balances  pressure-gradient  acceleration  in 

large-scale  flow;  mass  field  determines  velocity  field 
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Hydrostatic  Balance:  gravity  balances  vertical  pressure-gradient  in  large-scale 

flow;  temperature  determines  mass  field 
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Atmospheric  Governing  Equations 


■  Conservation  of  Mass:  Continuity  Equation 

^  =  -V-(pV) 

ct 

■  Conservation  of  Energy:  Thermodynamic  Equation 
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Atmospheric  Governing  Equations 


■  Equation  of  State 

p  =  pRT 

■  Conservation  of  water  vapor  mixing  ratio 


^^  =  -V(pyq)+p{E-C) 

ct 


Total  amount  of  water  vapor  in  a  parcel  Is  conserved  as  the  parcel  moves  except 
where  there  are  sources  (evaporation  E)  and  sinks  (condensation  C) 
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Equations  Summary 


■  Seven  Equations  with  seven  unknowns 
■  u,  V,  w,  T,  p,  p,  q 
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Introduction  to  Meteorological 

Observations 

■  In  situ 

■  Surface 

■  Upper-air 

■  Buoys 

■  Aircraft 

■  Ships 

■  Remotely  sensed 

■  Radar 

■  Wind  Profiler 

■  Satellite  soundings 

■  Satellite  winds 

■  May  not  be  meteorological  variables  (i.e.  radiance, 
reflectivity,  etc) 
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Global  Surface  Land  Observation  Coverage 


Global  Surface  Ocean  Observation  Coverage 


Global  Radiosonde  Coverage 


Global  Aircraft  Coverage  -Typical  6  hour  period 
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Global  Satellite  Wind  Coverage 


Polar  Orbiting  Sateiiite  Temperature 

Sounding  Coverage 
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Introduction  to  Numerical  Weather 

Prediction 


■  Global  models-  Spectral 

■  Spherical  harmonic  representation  of  data 

■  Equations  of  motion  in  amplitude  space 


■  ECMWF-  European  Center  for  Medium-Range  Weather  Forecasts 

■  T799L91-  Horizontal  grid  length  of  25km  and  91  vertical  levels 


■  GFS- Run  by  NCEP  (NOAA) 

■  T382L64-  Equivalent  to  about  40km  resolution  with  64  vertical  levels 


■  NOGAPS-  Fleet  Numerical  (Navy’s  Global  Model) 

■  T239L30-  55km  with  30  vertical  levels 


Friday  15  JUne  2007  OOUTC  OECMWF  Forecast  U072  VT;  Monday  18  June  2007  OOUTC 
Surface:  Mean  sea  level  pressure  /  850-hPa  wind  speed 
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Introduction  to  Numerical  Weather 

Prediction 


■  Limited  Area  (Regional)  Models-  Finite  Difference 

■  Nested  within  larger  global  or  regional  models 

■  Boundary  conditions  via  global  model  or  bigger  regional  model 

■  WRF- As  run  by  NCEP  (NOAA) 

■  12km  resolution  at  60  vertical  levels 

■  MM5-  As  run  by  AFWA 

■  Nested  Grid  45km,  15km  resolution  at  42  vertical  levels 

■  5km  output  over  selected  regions 


MM5/WRF  GLOBAL  WINDOW  CONFIGURATIONS 


LAST  MODIFIED;  Thu  Jun  14  18:15:31  2007 


11-JUI-07  https://weather.afwa.af.mil/data_links/MWORLDLOCA.GIF 


NWP  Forecast  Skill  Improvement 


100 


Year 


Anomaly  Correlation  ofSOO-hPa  height  for  3-,  5-,  and  7-day  forecasts  for  the  ECMWF 
operational  model  as  a  function  of  year.  Top  and  bottom  of  each  band  correspond  to  Northern  and 

Southern  Hemispheres,  respectively. 


■  Difference  between  hemispheres  has  neariy 
disappeared  in  the  past  few  years  due  to  the 
successful  assimilation  of  satellite  data 
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50hPa  —  lOOhPa  -  200  hPa  -  500  hPa  - 1000  hPa 


Anomaly  Correlation  Coefficient  for  ECMWF  forecasts  for  different  levels  over  the  Northern 

Hemisphere  in  2004. 
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Introduction  to  Numerical  Weather 

Prediction 


■  Sensitivity  to  parameterization  of  sub-grid  processes 

■  Introduction  of  chaos  and  error  into  models 

■  Model  has  0(10^)  values  to  calculate  (degrees  of 
freedom)  and  observational  quantities  are  0(10^) 

■  Problem  statement:  NWP  is  an  under-determined 
initiai  vaiue  probiem;  how  do  we  obtain  the  most 
reaiistic  representation  of  the  initiai  condition? 
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Three  Atmospheres 


We  deal  with  three  atmospheres 


■  Real  atmosphere 

■  Unknown 

■  Observed  atmosphere 

■  Data  coverage  gaps 

— vertical,  horizontal,  temporal 

■  Observation  error 

— random  and  systematic 

■  Analysis/model  atmosphere 

■  Observation  limitations 

■  Data  assimilation  system  limitations 

■  Model  limitations 
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Three  Sources  of  Information 


IMPERFECT 

OBSERVATIONS 


IMPERFECT 

FORECAST 

BACKGROUND 


ERROR  STATISTICS 
DESCRIBING 
IMPERFECTIONS 


BEST  LINEAR  UNBIASED 
ESTIMATE  OF  ATMOSPHERIC 
STATE 


Analysis/Forecast  Cycle 


The  forecast  model 
provides  the 
background  estimate 
of  the  current 
atmospheric  state 


The  analysis  provides 
the  initiai  conditions  for 
the  next  forecast 


These  corrections  are 
added  to  the 
background  to  form 
the  anaiysis 
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The  data  assimilation 
algorithm  produces 
corrections  to  the 
model  fields 
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Methods  of  Data  Assimilation 

■  Optimal  Interpolation  (01) 

■  3-Dimensional  Variational  Assimilation  (3DVAR) 

■  4-Dimensional  Variational  Assimilation  (4DVAR) 

■  Extended  Kalman  Filter 

■  Ensemble  Kalman  Filter 
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Methods  of  Data  Assimilation 


■  Optimal  Interpolation 


number  of  model  number  of 
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Methods  of  Data  Assimilation 


■  Optimal  Interpolation 

■  Optimal  analysis  is  found  by  minimizing  the  analysis  error 
variance  (A)  by  finding  the  optimal  weights  of  observation 
increments  through  a  least  squares  approach 

=  analysis  error  variance 
Gq^  =  observational  error  variance 
G^^  =  background  error  variance 
W  =  optimal  weight 
W  =  weight  matrix 

B  =  background  error  covariance  matrix 
R  =  observations  error  covariance  matrix 
H  =  linear  observation  operator  matrix 
I  =  identity  matrix 

A  =  analysis  error  covariance  matrix 


Minimized 


1  1 
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A  =  (I-WH)B 
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Methods  of  Data  Assimilation 


■  3DVAR  Introduction-  Cost  Function 

Why  use  the  term  Cost  Function? 


Xg  =  analyzed  value 
y  =  observed  value 
=  background  value 

Think  of  the  function  as  defining  how  far  away  the 
anaiysis  vaiue  is  away  from  the  input  vaiues.  You  wiii  be 
“charged”  for  not  fitting  the  observation.  You  wiii  aiso  be 
“charged”  for  not  fitting  the  background.  You  choose  the 
anaiysis  to  minimize  your  “cost.” 
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Methods  of  Data  Assimilation 


■  3DVAR  Introduction-  Vector  Form  of  Cost  Function 

Vector  Form: 

Minimizing  J(x)  with  respect  to  yields: 

^a=^b+  + R)”'  (y-^b) 
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Methods  of  Data  Assimilation 


■  3DVAR 

■  Minimizes  cost  function 

■  Performs  anaiysis  at  set  times 

■  Only  observational  data  from  set  times  are  inciuded  or  observations  in  a 
+/-  time  window  are  lumped  together  and  given  essentially  an  equivalent 
time 

■  Employs  a  forward  operator  or  observation  operator:  H 

■  Interpolates  observations  spatiaily  to  grid  points 

■  Converts  observed  quantities  (i.e.  sateilite  radiances,  radar  refiectivities) 
into  modei  variabies  (i.e.  temperatures,  humidites) 

J(x)  =  + 

I  ' - V - '  ' - V - ' 

Distance  to  forecast  Distance  to  observations 

At  analysis  time 
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Methods  of  Data  Assimilation 


m  4DVAR 

■  Minimizes  cost  function 


■  Cost  function  is  weighted  difference  between  modei  forecasts  during 
an  assimiiation  window  and  coincident  observations 


J(x)  =  min— 

2 


N 


(xito  )  -  (to  ))B“‘  (Xj  (to )  -  x„  (to  ))  +  X  (>”/  -  H*/  (>'/  -  ) 


Distance  to  background 
at  initial  time 


i=0\. - ^ ^ - 

Distance  to  observations 
in  a  time  window  interval 


■  4D-Var  seeks  initiai  conditions  such  that  the  forecast  best  fits  the 
observations  within  the  assimiiation  intervai 


■  Computes  the  increments  at  observation  times  during  a  forward 
integration  using  the  forecast  modei  and  then  integrates  these 
weighted  increments  back  into  the  initiai  time  using  the  adjoint  modei 

(L^) 
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Methods  of  Data  Assimilation 


X 


assimilation  window 


4DVAR  data  assimilation  window 
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Methods  of  Data  Assimilation 


■  Extended  Kalman  Filter 

■  Same  weighted  approach  as  Oi,  but  using  a  background  error  covariance 
matrix  (P)  that  evolves  with  the  forecast  rather  than  a  constant  background 
covariance  matrix  (B) 

■  Forecast  error  covariance  is  obtained  using  a  tangent  linear  model  (L)  to 
transform  the  perturbation  from  the  initial  time  to  the  final  time  and  the 
adjoint  model  (L''')  to  “advance”  the  perturbation  backwards  from  the  final 


to  initial  time  to  optimize  the 

initial  conditions 

Forecast  step: 

n 

P„ 

=  L„  A„_,L^  +  Q„ 

Analysis  step: 

< 

=  xi  +^^,Xyn-Hxi) 

The  optimal  weight  matrix  : 

K, 

=  p,(r+hpX)' 

The  analysis  error  covariance: 

A„ 
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Methods  of  Data  Assimilation 


■  “Extended  Kalman  filter  is  gold  standard  of  data 
assimilation”  (Kalnay,  2006) 

■  A  poor  initial  guess  can  be  transitioned  through  time  to 
provide  the  best  linear  unbiased  estimate  of  the  state  of 
the  atmosphere  and  its  error  covariance  provided 
observations  are  frequent  and  system  is  stable  (Kalnay, 
2006) 

■  Very  computationally  expensive 

■  So  can  it  be  done  more  cheaply? 
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Methods  of  Data  Assimilation 


■  Ensemble  Kalman  Filter 

■  Estimates  forecast  error  covariance  matrix  from 
ensembie  of  forecasts  initiaiized  by  random 
perturbations  added  to  the  same  sets  of  observations 

■  Computationai  cost  increased  0(10^)  from  Oi  or  3DVAR, 
but  cheap  compared  to  Extended  Kaiman  Fiiter  whose 
cost  is  0(10^) 

■  Tangent  iinear  or  adjoint  modei  not  necessary 

■  Most  promising  approach  for  the  future 
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Conclusions 


■  Limit  of  predictability  is  ~2  weeks  due  to  sensitive 
dependence  on  initial  condition  (chaos) 

■  Quality  of  forecast  therefore  is  critically  dependent  on 
the  quality  of  the  initial  condition 

■  Throughout  50  years  of  NWP,  methods  of  obtaining 
that  initial  condition  have  evolved  with  increasing 
sophistication  of  NWP  models,  new  remotely-sensed 
observation  types,  and  increasing  computing  power 
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Questions? 
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